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I. INTRODUCTION 

The existence of gravitational waves (GW), which is predicted in the theory of general 
relativity, has long been verified 'indirectly' through the observations of Hulse and Taylor 
HI]. The inspiral of the members of the binary pulsar system named after them has been 
successfully accounted for in terms of the back reaction due to the radiated gravitational 
waves Il|,@]. However, detection of such waves with man-made 'antennas' has not been 
confirmed so far. Nevertheless, this problem has received a lot of attention this decade, 
especially, due to arrival of laser-interferometric detectors, which are touted to have the 
sensitivity required for detecting such waves. 

In the past, a sizable amount of research has been done on the problem of detecting 
gravitational waves using a single bar or interferometric detector. However, very little work 
has been devoted on developing techniques to optimally analyze the data from a network 
of such detectors to seek the presence of coalescing binary signal. As has been argued in 
the past (see, eg., Ref 0), for a given false-alarm probability, the threshold for detection is 
lowered as the number of detectors is increased. This increases the probability of detection 
by a network rather than a single detector, provided the observer accepts only coincidences. 
One of the early papers which came close to discussing the problem of detection of these 
waves using a network was that of Finn and Chernoff [Q. This paper observed that since 
the orientations of the two LIGO detectors were very similar, their joint sensitivity was 
larger than any one of them. Another work which dealt with the issue of detection using a 
network was that of Bhawal and Dhurandhar 0. The main aim of this paper was to find the 
optimal recycling mode of operation of the planned laser interferometric detectors for which 
a meaningful coincidence detection of broadband signals could be performed. However, 
none of these earlier papers addressed the issue of how a network of detectors with arbitrary 
orientations can be optimally used as a "single" detector of sensitivity higher than that of 
any of its subsets of individual detectors. One of our main aims is to show precisely how this 
can be achieved. In the process, we will arrive at a network statistic based on the individual 
detector outputs that can be used to ascertain the presence of a signal in them with a given 
level of confidence. 

We note that the use of a network has nevertheless received considerable attention in 
the context of the parameter estimation problem. Some of the notable works that address 
this issue are Refs. []6|-|10|. The prime motivation in the use of networks in this regard is 



that the larger the number of detectors, smaller the errors in estimated values of the binary 
parameters. However, the starting point in these approaches is the assumption that the 
problem of detection has already been addressed and the detector specific chirp filters that 
result in "super-threshold" cross-correlations with the individual detector outputs, have been 
picked. 

Here, we formulate the data analysis problem in the case of the coalescing binary signal 
for a network of, say, A^ number of laser interferometric gravitational wave detectors that 
have arbitrary orientations, but are located at the same site. The noise in each detector is 
assumed to be additive and Gaussian. Also, the noises in different detectors are taken to 
be independent of one another. We use the maximum likelihood method for optimizing the 
detection problem. 

The paper is organized as follows. In Sec. |I|, we define the various coordinate frames, 



such as the detector frame and the wave frame, that we use in our calculations. We describe 
some known representations of the Newtonian signal corresponding to gravitational waves 
from a coalescing binary. In Sec. |IT1|, we present a new representation for this signal in 
terms of the complex expansion coefficients of the wave and the detector tensor in a basis 
of STF tensors of rank 2. Section |V| shows how the detection problem can be optimally 
addressed using the maximum likelihood method. In Sec. 0, we present the analysis for 
the improvement of the sensitivity of a network as a function of N. Finally, in Sec. [V^, 
we discuss how our results continue to hold for the restricted post-Newtonian waveforms. 
We also mention how our formalism can be extended to address the detection problem for 
a network of spatially well separated detectors. 

We use the following convention for symbols in this paper. Variables characterizing 
the network are displayed in the Sans Serif font. A parenthetic index in the superscript or 
subscript of a variable identifies a particular detector. Network- or individual detector-based 
variables that are complex are denoted by uppercase letters, whereas the lower case letters 
are reserved for real variables. Note that quantities such as the gravitational constant, G, 
though written in upper case, are not complex since they do not represent any inherent 
characteristic of the network or an individual detector. Also, we define the complex inner 
product as {AA(^j^, Bf) = A*B{A(^j^, f). By our convention, all the quantities featuring 
in the above expression, except f, are complex. Moreover, A(7) denotes a variable that 
characterizes the J-th detector in the network, where / is a natural number. Also, f denotes 
a network-based variable. 



II. PRELIMINARIES 

We describe the various coordinate frames in terms of which we will analyze the different 
polarizations of an incoming wave. Let (X, Y, Z) be the orthogonal Cartesian coordinates 
connected with a weak plane gravitational wave traveling along the positive Z-direction; X 
and Y denote the axes of the polarization ellipse of the wave. Let (x, y, z) form a right- 
handed coordinate system that describes a fiducial detector (henceforth referred to as the 
"fide" or the "network frame" ) . Let us define the Euler angles 9 and to give the incoming 
direction of the wave, and ip to denote the angle between one semi-axis of the ellipse of 
polarization and the node direction. The orthogonal matrix transformation from the wave 
frame to the fide is thus defined by the Euler angles {(f), 9, ip}. The orthogonal matrix 
transformation from the fide to the frame of the J-th detector is defined by the Euler angles 

{«(/), A^)'7(/)} • 

A gravitational wave is represented by metric tensor fiuctuation, hij, about the vacuum. 
In the transverse trace-free gauge, its non-vanishing components in the wave-frame are 
hxx = —hyy=h+, hxy = hyx=hx). Here, h^ and hx are the two polarizations of the waveform. 
In the Newtonian approximation, they are: 



2A/'a-i/^(t) 1 + cos- e 



/i^(t) = niZ ^ ^ ' 7' cos[x(t) + ^] , (2.1a) 



/i^(t) = ±^li^ \LLcosesm[xit)+6]. (2.1b) 



Above, A/" = 2G^/'^M.^f'^{'KfaY^'^/c^ , r is the luminosity distance from the earth to the 

binary, J\4. is the "chirp" mass defined by A^ = (1 + z)iJt'/^m?^^ ^ where m = mi + m2 is the 
total mass of the binary, /i is the reduced mass, z is the cosmological redshift of the binary, 
and c is the speed of light in vacuum. The angle e is the angle of inclination of the binary, 
i.e., the angle between the line of sight and the vector normal to the orbit of the binary, and 
6 is an initial phase of the orbital motion. The frequency of the gravitational wave is twice 
the orbital frequency and is given by 



f{t;ta,M) 




3/8 



TV \256 M^/^tc{ta,M)-t] 



(2.2) 



where ta is the time of arrival of the signal (such that f{ta) = fa = lOHz) and tc is the time 
at which coalescence occurs. Inverting the above equation after setting f{ta]ta,Ai) = fa, 
we get the time of coalescence: 



tc{ta,M) =ta 



.3\5/3 



256A45/3(7r/J8/3 \^^^ 



(2.3) 



Finally, x{t) = 27i ^ f{t')dt' and a{t) = 1 - (t - Q/^, where ^ = 
3.00(A4/Mq)~ ' (/a/100 Hz)~ ' sec, is the chirp parameter. Note that a total of eight 
independent parameters, viz., {r,6,9,(f),ip,e,ta,^} are required to specify this signal. The 
ranges of the four angles are as follows: 6* G (0, vr), G (0, 27r), tp G (0, 27r), and e G (0, vr). 
It can be shown that the signal at the fiducial detector is []TT| 



i{t) = 2Ka{t)-^/^ cos (z/„e (l - a{0^) + 6 - r]) 



where 
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1 + cos^e 
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cos^ e 



tanr] 



2cose 
1 + cos^ e 



(2.4) 

(2.5a) 
(2.5b) 



with 7] G (— 7r/4, vr/4). The signal at the J-th detector can be expressed in terms of the 
quantities defined above as 



Hi) 



(^) = 0(i)+h(^i)+{t) + 0(^i)^h(^i)^{t) . 



(2.6) 



where /i(/)+ and /i(/)x are the two polarizations of the wave arriving at the J-th detector. Also, 
0(/)+ and 0(7) X are the beam pattern functions, which depend on {0, 6, ip} and the orientation 
angles {q;(/),/3(/),7(7)}. The problem with the above representation for the signal is that it 
mixes up the factors dependent on the detector specific Euler angles {«(/), /5(/), 7(/)} and 
those dependent on the angles {0, 9, ■?/'}. We now give a different representation of the signal 
where this problem does not occur. It will prove useful to address the detection problem for 
a network. 



Consider the complex null vector M = l/-\/2(ex + iGy), where e^ and ey are real unit 
vectors in the X and Y directions, respectively, of the wave-frame. Then the wave tensor 



Wij is defined as 



Wij = h+Re{MiMj) + h,,lm{MiMj) , (2.7) 



which is a real symmetric trace-free (STF) tensor. The components of M in the detector axis 
are 0: M = l/-\/2(cos (f) — i cos 6 sin 0, sin + z cos 6 cos cj), i sin 6) exp{—i^p). A detector can 
also be represented as an STF tensor. For an interferometer with arms along the directions 
111 and n2 (both being unit vectors), the detector tensor d is given by 

dij = riuriij - n2in2j . (2.8) 

The response amplitude of the detector or, equivalently, the signal is just the scalar product 
of the wave and detector STF tensors, 

s = w'^dij , (2.9) 

where it is implicit that the Einstein summation convention holds over the repeated upper 
and lower indices i and j. In the following analysis, we will specifically consider a network 
of interferometric detectors. However, the generalization to bar detectors is straightforward. 
Since we extensively deal with STF tensors of rank 2, enunciating some frequently used 
properties of such objects is in order. Q Since such tensors have five independent elements, 
they can be expanded in terms of (location-independent) "STF-/" tensors, ylli, with rank 
/ = 2. They are related to the spherical harmonics as follows: 

Y2rn{0,^) = y2'mninj, where m = ±2, ±1, , (2.10) 

where n = (cos sin 6', sin sin 6*, cos 6). There are five independent 3x3 complex matrices, 
ylln, obeying the normalization 

yLyfr'* = ^^:^ ■ (2-11) 

In this paper, we will also be interested in the behavior of these STF tensors under action 
of an element g{a,P,'y) of the rotation group S0(3), where (a,/5,7) are the Euler angles. 
Consequently, we mention that under such an action, the spherical harmonics obey the 
following transformation law: 

Y2m'{e',<P')= Yl T_^-'^{a,f3,^)Y2Ud,<P) , (2.12) 

m=-2 



where T_m' "^ are the Gel'fand functions of rank 2 |12,B|. 



^For a detailed discussion, see Refs. [12,13|. For a more selective reading for immediate use, we 
refer to Ref. H. 



As shown in Ref . [^] , in the detector frame the wave tensor Wij can be expanded in terms 
of the STF-2 tensors as: 




w'^{t) = W— [/,+ (t)(T2" + T_2") - ^/^x(t)(T2" - T_2")] yi , (2.13) 



where the expansion coefficients are combinations of the Gel'fand functions, which depend 
on the parameters (0, 6, ip). For interferometric detectors with arms making an angle of 2^2, 
the only non- vanishing detector-tensor components in its own frame are di2 = d2i = sin2f2. 
The following analysis, where we will deal with the case Q = 7r/4, can be easily generalized 
to other values of Q. 

Consider the J-th detector of a network. Using the wave- and detector-tensor components 
in Eq. ( |2.9|) , it was shown in Ref. that the signal takes the form: 

S(,)(t) =J(7)„T„^(0,^, V^)r(,)r («(,),/?(,), 7(/)) , (2.14) 

Note that T(/)m"' is to be distinguished from T^^ in that T^" = Tm"'{(f),0,ilj), whereas 
7'(/)m" = T„"(a(/),/3(/),7(/)). Above, 

/(7)„ = ±^(/i+(,)-2/ix(/)) = ±J(/) for m = 2, n = ±2 , (2.15a) 

/(7)„ = ±^(/i+(/) + ih^n)) = TJiD for m = -2 , n = ±2 . (2.15b) 



The signal given in ( p.l4|) has the advantage of keeping factors dependent on the two sets 



of Euler angles separate. 



III. A NEW REPRESENTATION FOR THE SIGNAL 

To address the detection problem, we will be directly dealing with a construct known 
as the likelihood ratio (LR) (to be defined in Sec. fV\j . The LR is a non-linear functional 
of the signal. To keep this functional form simple we develop a new representation for the 
signal based on ( p.l4|) . 

In terms of J(/), the signal given in Eq. ( |2.14| ) takes the form 



S(,) = 2Re(zJ(*,)r(,)) , (3.1) 

where !(/) = i (t^j)^ - T^j^l^ T.s". Next, we define: 

hcii) = a'^j]^^{t)cosx{i)it) = g{i)S(i)o , 

hs{i) = a'^j]^^{t)smx(i){t) = 5'(/)-S(/),r/2 • (3.2) 

where ^'(7) is the maximum signal-to-noise ratio for the J-th detector obtainable by using an 
optimal filter: 

2 _ /" rj. \'''C,s[j)\ /„ „N 

9(1) = J ^f^iiy^ ' (3.3) 



s|^ (/) being the noise power-spectral-density (p.s.d.) of the J-th detector. Above, S(/)o and 
S{i)tt/2 can be represented as the real and imaginary parts of a complex quantity: 

S{i){t) = S(/)o + is(i)^/2 , (3.4) 

which has a norm equal to 2. 

Let us now express r(7) = 7(7)0 -|- «7(/)7r/2, where 7(7)0 and 7(7)77/2 are, respectively, the 
real and imaginary parts of r(7). Define 

W(^i) = 9{i) (7(7)0 cos 7] - 27(7)V2 siu 77) = W(7)e"*'^(-') , (3.5) 

where the last expression is the polar form of W(^iy Note that the only signal parameters on 
which 14^(7) depends are {(j),9,ip,v}- Armed with these definitions, the signal in Eq. ( p.l[ ) 
can be re-expressed as 



S(7)(t)=2«:Re(l^(*,)%)) , (3.6) 

where we have defined 

R(i)it) = r(7)o + ir (1)^/2 = S{^i)e'^ , (3.7) 

which is just the rotated S'(7). Here r(7)o and ^(7)77/2 are, respectively, the real and imaginary 
parts of -R(7). Equation ( p.6|) is a new representation for the signal that we will find useful 
in obtaining the maximum likelihood ratio below. 

We end this section by arriving at a relation between the complex variable W(i) and the 
detector tensor. First, by using Eqs. ( p. 131) and ( p.7|) we obtain the inner products between 



the wave tensor w^^ and the real and imaginary components of -R(7): 

2«:^(7)a*^ = (w*^(t),r(7)o) and 2«;^(7)/5*^' = (u7^^(t), r(7)V2) (3.8) 

where we have defined two new STF tensors aij and f3ij. These are real functions of 

{0,6',V',r7}. 



By comparing the different representations (|2.9| ) and (|3.6| ) of the signal we obtain 



W(i) = A'^ (^(7)t/(7)i,) , where A'^ = {a'' + t(3'') , (3.9) 

A*-^ is a complex STF tensor dependent on {(p, 6,1^,1]}, when expressed in the fide frame. 
Note that 

|iy(7)r = ley'' {9(i)d(i)r,) r + \P'' {9ii)d(i)r,) r ■ (3.10) 

Above, up to an r-dependent factor, |Vr(7)P can be interpreted as the total power transferred 
to the J-th detector. More appropriately, it is the gain factor associated with the /-th 
detector. We have resolved it as a sum of the fractions of power transferred to the detector 
by the two polarizations, respectively. It can be shown that, up to a factor of (^(7), W(^j) is 
just a direction cosine that is dependent on the set of angles {(f),9,ip,r]}. This is what one 
would expect from the above interpretation of W(i) as a gain factor. 



IV. ADDRESSING THE DETECTION PROBLEM FOR A NETWORK 

The signal from a coalescing binary will typically not stand above the broadband noise of 
the interferometric detectors; the concept of an absolutely certain detection does not exist in 
such a case. Only probabilities can be assigned to the presence of an expected signal. In the 
absence of prior probabilities, such a situation demands a decision strategy that maximizes 
the detection probability for a given false alarm probability. This is termed as the Neyman- 
Pearson criterion (see, eg., Ref. fl^)- Such a criterion implies that the decision must be 
based on the value of a statistic called the likelihood ratio (LR). It is defined as the ratio of 
the probability that a signal is present in an observation to the probability that it is not. 

For a network of detectors we obtain this statistic as follows. We assume that the 
noise at each detector is additive, Gaussian, and both statistically as well as algebraically 
independent of the noise in any other detector in the network. Under these conditions, the 
network LR, denoted by A, is just a product of the individual detector LR's. Similarly, the 
logarithmic likelihood ratio (LLR), In A, can be verified to have the same form as for an 
individual detector, namely, [p^pj. 

InA = (s,x)ArH^ - -(S, s)ArH/ , (4.1) 

where the normalized set of signals are denoted by a single network vector 

sW = (s(i)(t),S(2)(t), ,S(^)(t)) , (4.2) 

A^ being the number of detectors in the network. The subscript NW denotes that the inner 
product is defined on the network space. Similarly, the individual detector outputs X(/)(t) 
are combined to form the network vector 

x(t) = (^X(i)(t),a;(2)(t), ,X(^N){t)) . (4.3) 

Thus, in terms of the individual detector signals, the LLR is 

N 2 TV 

InA = ^(S(7),X(/))(7) - - I](S(/),S(/))(7) 
7=1 -^ 7=1 

TV ^ 

= bY.{Hi),Xii))ii)~-b' , (4.4) 



where 



7=1 



b = 2fi; 



\ 



N 

E«^(7) , (4.5) 

7=1 



is the norm of s and Z(/) = S(/)/b. The aim now is to maximize the LLR over all eight 
parameters to obtain the maximum (logarithmic) likelihood ratio (MLR). It is the MLR 
that must be compared with a threshold value to ascertain the presence or absence of signal 
in the detector output, with a given level of confidence. 
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We now analytically maximize the above expression with respect to as many of the eight 
parameters as possible. Note that the luminosity distance r appears only through b in LLR. 
Maximizing it with respect to b yields 

2 ^ /AT \ 2 

(4.6) 



N 



2 
where we have defined 



lnA|b = - (E(2;(/),X(,))(,)) =- 



v/=i 



Re(e-^^EQ(,)q,)) 



Wm 
<3(/) = ^^^=^= and C(*,) = (S(7),X(7))(7). (4.7) 

VS/=i^(/) 

The network vector S, with the >S'(/)'s as its components, is the matched network- filter. 

Next we maximize the LLR in Eq. (|4.6| ) with respect to 6 for, apart from the phase 
factor, none of the other terms there depend on it. This gives 

2 



InAk 



b,5 



2 



N 

(4.8) 



1=1 



which is a function of six parameters, namely, {(f),6,tlj,ri,ta,^}- Note that when all the 
detectors are "closely" located, it is only the Q(/)'s that depend on four angles {(f),9,'ip,ri}] 
the C(/)'s then depend only on {ta,C,}, with all the times of arrival being equal. We will 
refer to this situation as the "same-site" approximation. When the detectors are spatially 
well separated, the C(7)'s will depend on {4>,9} as well. 

To obtain the MLR, we need to maximize over these remaining parameters. At this stage 
it is useful to define the surrogate statistic (SS), A' = InAj^^. For a network comprising of 
a total of A^ detectors located within a fraction of a wavelength, the SS is maximum when 
Q oc C, where C is the network correlation vector with C(/)'s as its components. Therefore, 
once C is known, the maximization procedure determines Q through the above condition 
and the fact that Q has a unit norm. However, the Q so determined will, in general, yield an 
overdeterministic set of equations for the four parameters {</>, 6, ip, r]}. On the other hand, if 
this set of equations can be solved to yield a physically realizable solution for the parameters, 
then the maximized LLR will have a simpler form: 

1 ^ A 

^ 1=1 ^ 

where C(/) is the magnitude of C(/). Above, A is a function of two parameters, namely, 
{ta,^}- Although A is a real quantity, we follow the established convention in literature to 
denote the LLR by an uppercase letter! 

It can be shown that the condition Q oc C is always realised for two detectors. Numerical 



calculations suggest that this result holds for three detectors as well [T^. However, for 
networks with a larger number of detectors we numerically find that this condition is not 
always realisable and one is forced to maximize A'l^^^, as given in Eq. ( [4.8[ ), over the four 
angles. Thus, for networks comprised of up to three detectors, the application of Eq. ( |4.9|) 
appears to be valid. We will limit our discussion to only such cases below. Hence, only the 
maximization of A over the two parameters, {ta,C,}, remains to be done. This is performed 
numerically along the lines of Sathyaprakash and Dhurandhar (see Ref. ||11|| ). 
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V. NETWORK SENSITIVITY 

To infer the presence of a signal from the outputs of the members of a network, one 
compares the value of the statistic A in Eq. ( |4.9| ) with a predetermined detection threshold 
Aq. As we show below, the value of Aq can be obtained (via the Neyman- Pearson decision 



criterion [|1J]) from the false alarm probability, Qo, associated with the event of detection of 
such a signal. For A < Aq, presence of a signal in the data is ruled out, whereas if A > Aq, 
then the detection of a signal in the data is announced. 

We now analyze the improvement in the sensitivity of a network over that of a single 
detector. Apart from the assumptions about detector noise mentioned in Sec. ^V\, we further 
assume that it is stationary, which implies that for the J-th detector we have 

{n(i){f)nUf)) = sm)6{f - f) , (5.1) 

where Sh{i) is the noise p.s.d. of that detector and the angular brackets imply ensemble 
average. In general, different detectors may have non-identical noise p.s.d. We will assume 
that the noise in any detector is white, i.e., Sh{j) does not vary with frequency. Note that it 
has zero mean, i.e, {n(i){f)) = 0, and its standard deviation is Sh{i)- 

Equation ( ^.91 ) shows that A is a sum of squares of independent random variables with 
Gaussian probability distribution functions (PDF). Thus, A itself must have the so-called 
X^ probability distribution. Hence, the PDF of A under the hypothesis that the signal is 
present. Hi, is: 

Pi(A) = 2(b^) exp[-(A + b2)/2]/^,_l(bv^) , (5.2) 

where A > and In-i{x) is the modified Bessel function of order (A^ — 1). Note that b^ 
is proportional to the total energy received from the source. When bv^ >> 1, the above 
expression approximates to 

PiiVX) = -^exp[-(v^ - b)V2] . (5.3) 

On the other hand, under the hypothesis that the signal is absent, Hq, the PDF of A is 

""f^' - '2(^-1)! ■ (^■^' 

which one can obtain from Eq. ( p.2| ) by taking the limit b — > 0. 
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FIGURES 




Threshold, A. 



FIG. 1. The plot of Qq as a function of the detection threshold Aq for different number, N, of 
closely located detectors in a network. 



We are now in a position to calculate the false-alarm probability: 



Qo 



Ao 



Po{A)dA = 1 - I{Ao/2VN, N - 1) 



where we have made use of the incomplete gamma-function 



I{u, p) 



u^/p+T 



x^e "^dx/pl 



(5.5) 



(5.6) 



For a given false-alarm probability Qo, Eq. ( ^.5|) allows us to obtain the detection threshold 
Ao of our statistic. Plots of Qq versus Aq for different values of A^ are given in Fig. 1. From 
these plots it can be inferred that Aq increases slowly with N. 

The detection probability Qa can be obtained by computing the area under the function 
Pi(yA) for VA > ^/A^. For bv^ >> 1, it is 



Qa 



p^iVA)diVA), 



(5.7) 



where pi is given by Eq. (|5.3| ). We now show that for a given Qd and Qo, the distance 
r up to which a network can probe increases with A^. This is tantamount to saying that 
the sensitivity of a network increases as a function of N. For simplicity, assume that the 
detectors are oriented in such a way that b^ is proportional to A^. This is the case when, 
e.g., the W(7)'s are all identical. Let Qd = 0.5, i.e., v^A^ = b. As A^ increases, Aq increases 
for a given Qq (see Fig. 1.). However, given the fact that b oc \/N/r, we have r oc ^/N/Aq. 
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TABLES 

TABLE L The ratio of the sensitivity of a network with A^ = 2, 3, respectively, relative to a 
single detector for Q^ = 0.5, and corresponding to different false-alarm probabilities. 



False-alarm probability, Qo 



No. of detectors, N 


0.33x10-^'^ 


0.67 xlO-^^ 


0.17 xlO"^^ 


0.33 X 10-^4 


2 


1.33 


1.33 


1.33 


1.34 


3 


1.55 


1.56 


1.57 


1.58 



Thus, given a specific binary, the ratio of sensitivity of a network to that of a single detector 
behaves as 



r{N) 



r{N = 1) 



A^ 



^ Ao(iV=l) ' 
. Ao(Ar) , 



1/2 



(5.8) 



which can be computed by using Ao(A^) from Fig. 1. These ratios, which are presented in 
Table 1. for A^ = 2, 3, clearly show that the sensitivity of such a network increases a little 
slower then yN . This implies an increase in the survey volume accessible to a network, 
which, in turn, implies an increased event rate. 



VI. CONCLUSION 

The problem of detecting inspiral waveforms from coalescing binaries via pattern- 
matching can be made more accurate by including post-Newtonian corrections in the corre- 
sponding filters. In this regard, it has been shown [|l^ that it is both necessary and sufficient 
to work with the restricted post-Newtonian chirp. The description of the resulting wave- 
form involves an extra parameter, apart from the set of eight parameters described above. 
However, as was shown by Sathyaprakash |T^ , for the astrophysically relevant range of pa- 
rameters, the effective dimensionality of the parameter space remains unchanged. Hence, 
even after the inclusion of the restricted post-Newtonian corrections in the filters, the de- 



tection problem can be addressed in the same manner as described in Sec. |V. 

When the detectors are spatially well separated, the cross-correlations, C(/)'s, will be 
dependent on the times of arrival, ta(/)'s, which are different from one another. Since a 
specific ta(i) depends on the location angles {0,^}, so will C(/). Hence the maximization of 
the SS over the four angles {0, 9, iJj, rj} that was performed for the same-site approximation 
above, can no longer be implemented in the present case. It can be shown that the SS can 
be recast in such a way that its dependence on the complementary angles {ip, t]} is isolated. 
This aids in the analytic maximization of the SS over these two angles. The maximization 
over the remaining four parameters {0, 6, ta, ^} can then be effected numerically on a four- 
dimensional parameter-space grid. Details of these calculations will be presented elsewhere 
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